Stable Skyrmions in spinor condensates 
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Globally symmetric spinor condensates in free space are argued not to support stable topological 
defects in either two or three dimensions. In the latter case, however, we show that a topological 
Skyrmion can be stabilized by forcing it to adopt certain density profiles. A sufficient condition 
for the existence of Skyrmion solutions in three dimensions is formulated and illustrated in simple 
examples. Our results pertain to Bose-Einstein condensation in 87 Rb. 



Topologically non-trivial configurations play a crucial 
role in our present understanding of ordering, dynamics, 
and criticality in basic models of statistical mechanics 
lj. Probably the best understood example is that of vor- 
tices and vortex-loops in the Ginzburg-Landau, or Higgs, 
0(2)-symmetric theory for a complex order parameter in 
two (2D) and three dimensions (3D), which are responsi- 
ble for the veryexistence of the disordered phase at high 
temperatures [2|- The next in order of increasing com- 
plexity is the 0(3)-symmctric Heiscnbcrg model, which 
allows texture-like Skyrmions in 2D y], and point- like 
hedgehogs in 3D 0. Whereas Skyrmions in 2D are un- 
derstood to provide an additional source of disorder for 
the low-temperature phase, the exact role of hedgehogs in 
the 3D Heisenberg model is already less clear 5J. Both, 
nevertheless, represent non-trivial spatial configurations, 
which by virtue of their distinct topology are locally sta- 
ble and separated from the ground state by an infinite 
energy barrier. 

In this Letter we study the issue of topological de- 
fects in the next simplest case: the order parameter with 
global 0(4) symmetry. Such a symmetry arises in com- 
plex condensates with an internal spin-l/2-like quantum 
number [(j, for example. Realizations of such spinor 
condensates are found in models of inflatory cosmology 
0, Bose-Einstein condensation of 87 Rb (8J, bosonic fer- 
romagnetism |g, Q, and in effective theories of high- 
temperature superconductivity [l(| . [Tj, and of decon- 
fined criticality Jl2j . The Higgs sector of the Weinberg- 
Salam model of electro-weak interactions represents an- 
other closely related example, with a spinor condensate 
coupled to gauge fields. It is easy to show that the trivi- 
ality of the first and the second homotopy groups of the 
three-dimensional unit sphere, S3, on which the order 
parameter lives at low temperatures, implies that there 
are no stable topological defects in two dimensions. 2D 
spinor condensates may therefore be expected to be in 
the disordered phase at all finite temperatures, precisely 
as described by the 0(4) non-linear tr-model (NLaM). 
In 3D, on the other hand, standard topological consid- 
erations suggest the possibility of a Skyrmion texture, 
thanks to the third homotopy group of S3 being the group 
of integers. Here, however, unexpected subtleties arise, 
which together with the growing relevance of the problem 
provide the motivation for the present work. 



It has been known that finite-energy textures in higher 
than two dimensions are generally unstable with respect 
to shrinkage |13J. Derrick's scaling argument, does not, 
however, forbid topological defects with infinite energy, 
hedgehogs in the 3D Heisenberg model being a prime ex- 
ample. The search for stable topological defects in the 
0(4) NLaM in 3D, however, has so far led to a neg- 
ative result [3, Ea, 0, [l3 • We first demonstrate 
that even with amplitude variations included the spher- 
ically symmetric Skyrmion remains unstable. This con- 
clusion follows from a useful mechanical interpretation of 
Skyrmion's differential equation, in which the radial de- 
pendence of particle's density appears as a 'source of dis- 
sipation' for a fictitious classical particle, which prevents 
it from oscillating with a full amplitude. This hindrance, 
however, may be turned into an advantage by forcing the 
same radial dependence to have a form that effectively 
serves as an 'energy pump'. This idea yields an integral 
condition on the density profile of a Skyrmion, which 
we illustrate with three qualitatively different analytical 
solutions of the non-linear Euler-Lagrange's equations. 
Skyrmion solutions are found, for example, for specific 
forms of the confining potential, which may be used in 
creating such configurations in the laboratory |g|. 

Let us begin by defining the system of two-component 
bosons in the continuum with the standard action in 
terms of the complex coherent states written as 

r? u 

S= / drd 3 f[$ t (9 r -V 2 -Ai+^(f))$+-($ t $) 2 ], (1) 
Jo 2 

where $t = ($*(:?, r), ^(z, r ))- We have also included a 
confining potential V(x), assumed a contact repulsion U 
for simplicity, and set h — 1 , and the boson mass 2m = 1 . 
[3 = 1/ksT is the inverse temperature. 

We are interested in finding classical, that is r- 
independent, locally stable field configurations of the ac- 
tion (1). The field $(a?) may be written in terms of its 
real amplitude f(x) and the normalized complex spinor 
a[x) as & = jq)^ w ith a^a — 1. The low-energy physics 
is then described by fluctuations of the spinor a, which 
can be identified with a four-component real vector of 
unit length. The symmetry of the order parameter is 
thus 0(4). Varying the above action leads to the requi- 
site Euler-Lagrange's equations 

f{a)V 2 a + (V 2 at)a)a = 2/V 2 a + 4(V/) • (Va), (2) 



V 2 /+[(Va T )-(Va)-l + F]/ + / d 



0. 



(3) 



where we have rescaled the lengths as x^/Ji — > x, the 
amplitude as fy/U/fi — > /, and the external potential 
as V/fi — > V. Solution needs to satisfy the boundary 
conditions 



lim j;7 z Va = 0, 



(4) 



which guarantee stability with respect to small rotations 
of the spinor a at the infinitely remote boundary of the 
system. 

Eqs. (2) and (3) form a set of five differential equations 
for real and imaginary parts of the spinor a and the am- 
plitude /. In a spherically symmetric potential V(x) we 
can partially solve them, however, by assuming the most 
general ansatz with the same symmetry J4j: f(x) = ,/(r), 

a'(x) = (sinw(r) cos9 + icosuj(r), sinw(r) sin^e 1 ^), (5) 

where x = (rsin#cos0, r sin sin 0, rcos9). Somewhat 
tedious but otherwise straightforward algebra shows that 
all four of Eqs. (2) will then be satisfied provided that 
the function w(r) satisfies the differential equation 



2 df(r) du{r) sin(2w(r)) 



d 2 uj{r) 2 
dr 2 r f(r) dr dr 

Similarly, Eq. (3) then reduces to 

d 2 f(r) 2df(r) , 



dr 2 



dr 



l<^> 2 + 



0. (6) 



(7) 



2( 



S -^l) 2 -l + V(r)}f(r) + f(r) = 0. 



The action for the ansatz (5) is 



S 

2(^f + V(r) 



dr 



((^) 2 
dr 



(8) 



.f 4 



]dr, 



so that Eqs. (6) and (7) may also be recognized as Eulcr- 
Lagrange's equations for the latter form of S. 

The ansatz for a(x) will be topologically non-trivial 
and wrap the S3, defined by the condition a^a = 1, N 
times if, as the radius varies from zero to infinity, the 
function u>(r) takes all the values from zero to Nir. Here- 
after we restrict our discussion only to the elementary 
Skyrmion with N = 1. 

If we turn the external potential off (V(r) = 0), it is 
easy to show that the Eqs. (6) and (7) do not actually 
have the Skyrmion as a solution, the assumptions to the 
contrary notwithstanding [l5j, |ly. To this purpose de- 
fine a new 'time' variable —00 < t = lnr < 00. Eq. (6) 
then becomes 



w(t) 



dW(uj) 
du) 



»?(t)w(t), 



(9) 
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FIG. 1: A numerical solution of Eqs. (6) and (7): w(r) 
(thicker curve) and f(r). 



with W(w) = (cos2w)/2, r\{t) = 1 + d\nf 2 (t)/dt, and 
d) = duj/dt. The function ui(t) may now be interpreted 
as a coordinate of a particle moving in the potential 
W(u>), coupled to a dissipative environment with the 
time- dependent coefficient of dissipation rj(t). Let us 
first neglect the amplitude variations dictated by Eq. 
(7). Dissipation then still persists, and r)(t) = 1. De- 
manding that w(0) = Eq. (6) implies w(r) ~ r for 
r < 1, so that the initial condition in our mechanical 
analogy is lo(— 00) = u>(— 00) = 0. This corresponds to 
rolling down from the top of the potential W(ui) without 
any initial kinetic energy. Independently of the slope of 
u>(r) at the origin the solution therefore always oscillates 
around and approaches w(oo) = 7r/2, dissipating energy 
and loosing its amplitude in the process. Eq. (6) for 
a fixed amplitude therefore admits only a meron (half- 
Skyrmion) as a solution. However, since when r ^$> 1, 
\uj(r) — tt/2| ~ 1/v^* f° r the meron, the boundary condi- 
tion in Eq. (4) is not satisfied. Consequently, the meron 
is unstable with respect to deformation u>(r) — > u>(r)—e(r) 
with e(0) = and < e(oo) < I 18]. 

Before relaxing the condition of fixed amplitude let us 
observe that solutions of Eq. (6) (still with f(r) = const) 
that satisfy the boundary conditions (4) and approach 
w(oo) = 7r as uj(r) ~ 7r ± 1/r 2 , do in fact exist. These 
solutions, however, are singular near the origin, w(r) ~ 
±l/r, and hence have an infinite action. We will refer 
to these as singular solutions and discuss their possible 
significance in a moment. 

Next, we retain V(r) = 0, but allow the amplitude 
variations. Finitcncss of the action, measured from the 
value for the ground state f(r) = 1, u>(r) — 0, requires 
then /(oo) = 1. Linearizing near the origin one still finds 
«(r) - Ar + 0(r 2 ), and /(r) = /(0) + (/(0)/6)(/ 2 (0) + 
3A 2 — l)r 2 + 0(r 3 ), where A and /(0) are finite con- 
stants. To have /(oo) = 1 we then find numerically that 
/(0) < 1, and f{r) to be a weakly and uniformly in- 
creasing function (Figure 1). The reason behind such a 
behavior is obvious: gradients of a(x) are largest near the 
origin, so suppressing the amplitude there somewhat al- 
ways lowers the energy. Energy density being integrable, 



however, the amplitude remains finite at the origin, and 
the Skyrmion is 'coreless'. With such a solution for f(r), 
however, one still finds w(oo) = ir/2 and Eq. (4) violated. 
This is evident from our mechanical analogy in Eq. (9) in 
which having the amplitude increasing with radius only 
adds to 'dissipation'. 

The stage is now set for the main part of the Let- 
ter. The preceding analysis points to a simple way 
to stabilize the Skyrmion. To have a solution of Eqs. 
(6) and (7) satisfying the requisite boundary conditions 
Lu(t = — oo) = d)(t = — oo) = 0, and u(t — oo) = w, 
u)(t = oo) = 0, it is necessary and sufficient that the total 
dissipation in our mechanical analogy vanishes. Mathe- 
matically [19J : 



r){t)uj 2 {t)dt = 0. 



(10) 



The total mechanical, and thus the potential, energy is 
then the same at the initial and the final 'time' t. The 
amplitude therefore must be a decreasing function at least 
for some radia, where it would 'pump' the energy back 
into the oscillator. This is our main result. In the re- 
maining we provide some specific amplitude profiles that 
indeed lead to stable Skyrmions. 

1) The trivial case is the one without dissipation: 
r)(t) = 0. This is equivalent to the particle density 



f(r) = -, 
r 

with A as a constant. This yields 



w(r) = 2cotf- 1 (r ) 



(11) 



(12) 



Inserting this solution into Eq. (7) gives the external 
potential that enforces such a density to be 



V(r) 



1 



A 



1 

4 r 2 



16r 2(V2-l) 

(1 + r 2 ^) 2 



(13) 



It is easy to check that this solution has a finite action, 
due to the singularity of the external potential at the ori- 
gin. The same singularity, however, makes this simplest 
example less than completely satisfying from a possible 
practical point of view. 

2) The difficulty in utilizing the condition in Eq. (10) 
is that the velocity <l>(t) of course depends on the dissi- 
pation r/(t) in a rather non-trivial way, so it appears that 
solving the differential equations (6) and (7) is always 
unavoidable in practice. This, fortunately, is not so, due 
to the following theorem: For potentials W(ui) that are 
even functions of lj — n/2, and dissipation coefficients 
r)(t) that are odd functions oft, solutions of Eq. (9) are 
such that 






(14) 



i.e. odd functions of time, measured from the moment 
when at the bottom of the potential. It is thus suffi- 
cient to have a density profile that implies an odd 77(f) to 
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FIG. 2: Topologically trivial ground state /o(r) (upper curve) 
and the Skyrmion's amplitude f(r), in the potential in Eq. 
(18) with B = l. 



stabilize the Skyrmion. As an illustration, consider 

n(t) =-tanh(t), (15) 

for which the the particle density is simply 

B 



f(r) 



1 



(16) 



with B as a constant. It is straightforward to check that 
the solution of Eq. (7) is then 



w(r) = 2cot~ 1 (r), 
and the corresponding external potential 

B 15 



V(r) 



1 



1 + r 2 (1 + r 2 ) 



2\2 ' 



(17) 



(18) 



Different density profiles can be similarly constructed, 
including those that would lead to a finite total number 
of particles. 

We may note that for the external potential in Eq. (18) 
there is yet another, topologically trivial, solution of Eqs. 
(6) and (7): w{r) — 0, fo{r). Numerical solution for /o(r) 
for B = 1 is plotted on Figure 2. Since fo{r) > f{r), the 
action, which at a stationary point takes the form 



4^a 



T-,=-\l r ' Z f»dr, 



(19) 



is actually lower for fo(r). This represents the ground 
state configuration in the potential (18). 

3) Finally, the Skyrmion may be stable in a constant, 
but discontinuous, particle density as well. Take f 2 {r) = 
e for r < 1, and f 2 (r) — 1 for r > 1, for example. Then 



d\nf 2 {t) 
~dt 



= -cs(t), 



(20) 



and the total dissipation will vanish for C tuned to 



C 



Mt) . 



-dt. 



(21) 



Such a discontinuity in the amplitude may be under- 
stood as providing a 'kick' in our mechanical terminol- 
ogy, which instantaneously increases the kinetic energy 
at t = so that there is just enough to reach the next 
top of the potential W(lu) at ui = it. This is possible be- 
cause for any regular solution of Eq. (6) (with / = const) 
which starts at Lu(r = 0) = there is a singular solution 
that intersects it at r = 1, and which asymptotically ap- 
proaches tt as the radius increases. Their difference in 
slopes at r = 1 determines the required injection of en- 
ergy measured by the discontinuity C. 

Previously, the Skyrmion was found to be stable nu- 
merically if the 0(A) symmetry of the action in Eq. (1) 
is broken down to 0(2) x 0(2), in a way that favors 
phase separation |17| . In contrast, our action is fully 
0(4) symmetric, and the stability obtains from the in- 
terplay between the form of the Skyrmion and its den- 
sity profile. Derrick's scaling argument is likewise evaded 
by the introduction of the preferred length scale for den- 
sity variation. We hope our results will help facilitate 
the production of these interesting topological objects in 
laboratory spinor condensates, such as a7 Rb. Recent ad- 
vances in (dynamic) manipulation of confining potentials 



seem particularly promising in this respect [2fJ. 

Finally, for completeness let us revisit the issue of de- 
fect's stability in 2D. Consider 

a f (x) = (sinAcos0(r) + icos A,sin Asin6>(r)e 10 ), (22) 

where x — r(coscj), sin(/>), and 9(0) = and 6(oo) = 
7r/2. For a fixed parameter A = n/2 this reduces to 
the configuration often discussed in literature 8]. It is 
evident, however, that by tuning A one may continuously 
deform the defect at A = n/2 to the trivial vacuum at 
A = 0, monotonically decreasing the action along the 
way. This construction exploits the fact that ^(S^) = 
1. Topological defects in spinor condensates are thus in 
general unstable in free space both in 2D and in 3D, 
although for rather different reasons. 
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